Ballistic aggregation for one-sided Brownian initial velocity 



Patrick Valageas 

Institut de Physique Theorique, CEA Saclay, 91191 Gif-sur-Yvette, France 



Abstract 

We study the one-dimensional ballistic aggregation process in the continuum limit for one-sided Brownian initial 
velocity (i.e. particles merge when they collide and move freely between collisions, and in the continuum limit the 
initial velocity on the right side is a Brownian motion that starts from the origin x - 0). We consider the cases where 
the left side is either at rest or empty at f = 0. We derive explicit expressions for the velocity distribution and the 
mean density and current profiles built by this out-of-equilibrium system. We find that on the right side the mean 
density remains constant whereas the mean current is uniform and grows linearly with time. All quantities show an 
exponential decay on the far left. We also obtain the properties of the leftmost cluster that travels towards the left. We 
find that in both cases relevant lengths and masses scale as and the evolution is self-similar. 

Key words: Adhesive dynamics, Ballistic aggregation, Inviscid Burgers equation, non-equilibrium statistical 
mechanics 



1. Introduction 

We consider in this article the continuum limit of a one-dimensional ballistic aggregation process, for the case 
of Brownian initial velocities (i.e. the initial velocity field is a Brownian motion). In such a model, point particles 
of identical mass m move on a line and perform completely inelastic collisions, that is, in binary collisions particles 
(or clumps) merge to form a single larger aggregate under conservation of mass and momentum (but dissipation of 
energy). Between collisions clumps move at constant velocity (free motion). Thus, without external forcing, the 
stochasticity is only due to the randomness of the initial velocities. This model was introduced in ijst], for the case of 
uncorrected initial velocities (i.e. white-noise case in the continuum limit), as a simple test-case for scaling arguments 
used in more general hydrodynamical or statistical systems. Indeed, this ballistic aggregation process can be seen as a 
simple model for the merger of coherent structures, such as vortices, thermal plumes, or cosmic dust into planetesimals 
within proto-planetary disks. 

In this context it is natural to investigate the late-time asymptotic scaling regime obtained for the case of uncor- 
rected initial velocities. Thus, one finds that the average cluster mass grows with time as fil^ with a large-mass tail 
for the universal mass distribution of the form e^™'/'" I^T^, 12]. When the number of particles is finite, at long times 



the system reaches a stationary "fan" state, where the velocities of the final clusters increase from left to right. This 
final state also shows many universal properties, such as the number and size distributions of final clusters and the 



size of the leftmost and rightmost clusters 11261 12111 . On the other hand, when the initial particle velocities are given 
by a Brownian motion this ballistic aggregation process can be related to a simple additive coalescent model (which 
does not take into account positions nor velocities), where each pair of clusters merges with a rate proportional to its 
total mass, independently of other pairs [4]. This also provides results for the statistics of dislocation of clusters in the 
time-reversed fragmentation process. 

In the continuum limit, m — > at fixed uniform initial density pQ, it is well known that this system can be mapped 
onto the Burgers equation in the inviscid limit 15, 23, 121. 



dv dv d^v ., 

^ + = f with e ^ 0+. (1) 

at ax ax^ 

Then, shock locations in the Eulerian velocity field v(x, f) describe particle aggregates of finite mass whereas regular 
points correspond to infinitesimal particles. It is clear that away from shocks Eq.O corresponds to free motion (in 
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the limit e — > 0^ where the right-hand side vanishes) and it can be shown that shocks conserve momentum 1^, which 
explains the relation with the ballistic aggregation process. 

The Burgers equation ([T]) itself is a nonlinear evolution equation that appears in many physical problems, such as 
turbulence studies I6j^20], the propagation of nonlinear acoustic waves [15], or the formation of large-scale structures 
in cosmology lfl4ll28ll . see the recent review |2] for a detailed discussion. In particular, the study of the statistical 
properties of the dynamics, starting with random Gaussian initial conditions, is also referred to as "decaying Burgers 
turbulence" in the hydrodynamical context |16], or as the "adhesion model" in the cosmological context |14]. In 
these frameworks, where the random initial velocity applies over all space, and may be homogeneous or only have 
homogeneous increments, one is interested in Eulerian quantities such as the velocity structure functions, the «-point 
velocity distributions, the matter density distribution, the mass function of shocks, or Lagrangian quantities such as the 
distribution of the displacement field. Then, it is customary to consider power-law initial energy spectra, Eo{k) oc k", 
where the initial velocity fluctuations scale as over size x (the white-noise case is « = and the Brownian case 

is n = -2). Then, for -3 < « < 1, a self-similar evolution develops [15, 23]. The integral scale of turbulence, which 
measures the typical distance between shocks and the correlation length, grows as L(f) ~ f2/(n+3) ^vhereas the tails of 

the cumulative shock distribution and velocity distribution satisfy ln[n(> m)] m"^^, ln[n(> |v|)] |v|"^-', for m — > 

cxj, |v'| — > oo, see [24, 22]. In such a context, the white-noise case, n - 0, corresponds to initial velocity fluctuations 
that are dominated by high wavenumbers, whereas they are governed by low wavenumbers in the Brownian case. 
Then, this latter case is of particular interest since in many hydrodynamical systems the power is generated by the 
larger scales. For instance, the Kolmogorov spectrum of turbulence, E(k) oc k^^^^, shows such an infrared divergence, 
whereas in cosmology the velocity fluctuations are also governed by scales that are larger than the nonlinear scales 
where structures have already formed in the density field. 

The connection with the Burgers equation ([1]) allows us to derive many results for the continuum limit of the 
ballistic aggregation process, taking advantage of its well-known Hopf-Cole solution ifisi lgll. (This also corresponds to 
the late-time evolution of the system if we keep a finite particle mass.) In particular, using the geometrical description 
of this solution in terms of first-contact points between the initial velocity potential and parabolas [6J, as recalled 
below in or the equivalent description in terms of the convex hull of the Lagrangian potential jQl, it is possible to 
derive closed analytical results for the specific cases of white-noise and Brownian initial velocities. Indeed, in these 
two cases the velocity or potential fields are Markovian which allows us to greatly simplify the analysis n, 271. 
Moreover, a specific property obtained for Brownian initial conditions is that the inverse Lagrangian map, x i-> q(x, t), 
where q is the initial Lagrangian position of the particle located at x at time t, keeps homogeneous increments at all 
times (on the right side for one-sided initial conditions and far from the origin for two-sided initial conditions) IsI IztIi . 
This can actually be extended to some Levy processes with no positive jumps ||3|l. 

Then, for both white-noise and Brownian one-sided cases, [5] found that the flux through the origin is a pure jump 
time-inhomogeneous Markov process and obtained its statistical distribution. While for the white-noise case clusters 
(shocks) are in finite number per unit length, which implies that the mass that has flowed to the the left side increases 
through a finite number of jumps per unit time, for the Brownian case clusters are dense (on the right side), which 
implies that on any time interval some small-mass clusters have crossed the origin [24, 5]. Some properties of the 
limit clusters travelling towards the left were obtained in lfl9ll for the white-noise case, while for the Brownian case it 
was found that shocks are dense to the right of the leftmost cluster, as on the right side, but no results for the statistics 
of the latter were obtained. Statistical properties of limit clusters (i.e. at the infinite time limit) were also obtained in 
II29II for initial velocities that are given by a Levy process with no positive jumps. 

For the white-noise case, [1121 \lM studied the late-time dynamics reached when the "excited" particles are re- 
stricted to the semi-infinite right side, or to a finite interval, and expand into empty space or a medium at rest. Many 
explicit analytical results can be derived in the continuum limit [13]. If the initial "excited" interval is finite and 
surrounded by empty space, the late-time evolution is ballistic and the characteristic length scales as L(f) oc t. Indeed, 
since the total mass is finite, the system eventually reaches a "fan" state with a finite number of clusters that move 
freely without anymore collisions. If the "excited" particles expand into a medium of uniform density at rest, the 
latter slows down the motion and the characteristic length scales as L(t) oc f'^^. For one-sided initial conditions, where 
the initial white-noise velocities apply to the semi-infinite right axis, one recovers the scaling law associated with the 
homogeneous turbulent case recalled above, L(t) oc t^^^ and M{t) oc fi/^. Thus, nontrivial mass density and current 
profiles develop over this scale, with a power-law tail p ~ |x|"^ on the far left. In addition, other nontrivial asymptotic 
mass profiles are obtained over scales La{t) oc t", with 2/3 < a < 1, that interpolate between the natural scaling oc fi/^ 
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and the ballistic regime cx t. This corresponds to forerunners that carry a mass oc f^^i"^^ interpolating between masses 
of order and of order unity. Finally, if the "excited" particles expand into a fiUed medium at rest, other profiles 
develop on the natural scale L(t) oc f^^^, but there is no more propagation on larger scales, a > 2/3, because of the 
slowing down by the left-side particles, and the density decays exponentially fast on the left side. 

In this article, we study how these results are modified when the initial velocities on the right side are given by 
a Brownian motion, instead of a white-noise spectrum, and we also consider the statistical properties of the leftmost 
("leader") cluster that has formed on the left side. As noticed above, the Brownian case is a template for large-scale 
forcing as opposed to small-scale forcing in the initial velocity field. Thus, we consider the case where the initial 
velocity field, vo(^) at time f = 0, is a Brownian motion on the semi-infinite axis q > 0, while the initial density po is 
constant over q >Q, and we write 

p(q,0)^po, Vo(?)= Cdq'^iq') and Mq)^ C dq' f dq" ^(q") owev q > 0. (2) 
Jo Jo Jo 

Here we introduced the velocity potential i//(x, t), with v = di^/dx, and a Gaussian white-noise ^(q), which we nor- 
malize by 

{^(q)) - 0, {^(q)^{q')) - D 6{q - q), whence (v'o(^)^) = Dq in the Brownian region, (3) 

where (..) is the average over all realizations of ^. In Eq.® we normalized the initial velocity and potential by 
vo(0) - and i/'o(0) = at the origin. We consider two cases, "F" and "£", where the left semi-infinite axis, <7 < 0, is 
either filled with particles at the same density po but with zero initial velocity (medium at rest), or empty (zero density, 
p = 0). Therefore, we complete the definition (|2]i by 

case "F" : p{q, 0) - po, vo(^) - and iAo(?) - over ^ < 0, (4) 
case"/?" : p{q,Q) - po, voiq) - and ifroiq) = v^q over ^ < 0, with v_ — > -oo. (5) 

Here, as in fl^, we used the fact that the empty case, "£", of (|5]l, can be obtained by keeping the same uniform initial 
density po over q <Q, while giving to these particles a velocity v_ that goes to -oo. Then, these particles immediately 
escape to the infinite left at f = 0^ and the Brownian particles with q > spread into empty space. These initial 
conditions can be summarized by stating that the initial velocity potential iJ/q is either continuous and constant out of 
the Brownian region (filled case "F") or goes to +oo (empty case "£"). 

At any point jc < on the left part, the system remains unchanged (at rest or empty) until some particles that 
originate from the right-side Brownian region have managed to travel down to position x. Note that once some 
particles have entered the left part they will keep travelling with a negative velocity forever However, their velocity 
can change as they may overtake slower particles or may be overtaken by faster particles that escaped at a later time 
from the Brownian region. Since particles do not cross, the leftmost cluster is associated with the particle, ^ = 0, that 
was initially at the left boundary of the Brownian domain. Once the latter has entered the left side, it keeps moving to 
the left and in case F it draws along all the matter that was initially at rest. 

In the context of the ballistic aggregation process studied in this article, and contrary to the hydrodynamical context 
where the Burgers equation ([TJ is used to investigate statistically homogeneous turbulence, the systems defined by 
Eqs.©-® are clearly statistically inhomogeneous and a current develops towards the left side as particles escape 
into the left part of the system and then keep travelling to the left forever. Therefore, we mainly focus on quantities 
that express this out-of-equilibrium propagation of matter towards the left. Indeed, the conditional probabilities to the 
right of the leftmost cluster are identical to the ones obtained for the two-sided Brownian-motion initial velocity ||27|] . 
In particular, the distribution of velocity increments and of the matter density are the same on the right side x > 0, see 
also H. 

We recall in section|2]the geometrical construction in terms of first-contact parabolas of the solution of Eq.([T]i and 
the associated Brownian propagators. Next, we first consider the case "F" of a filled left-side at rest, and we study 
the velocity distribution Px{v) as well as the probability p^hocked jj^^j. u^^^tg]- from the right side has already reached the 
position X on the left side by time t. We consider the mean density profile and current in section l3T2l Then, we derive 
the Lagrangian displacement field in section [33] and we obtain in section [341 the properties of the leftmost cluster. 
Finally, we consider the case "£" of the empty left side in section|4] 
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2. Geometrical construction and Brownian propagators with parabolic absorbing barrier 



As recalled above, in the continuum limit, m — > at fixed density po, the ballistic aggregation dynamics is fully 
described by the Burgers equation ([T]) in the limit of zero viscosity. As is well known |18l \%. substituting for the 
velocity potential t), with v = dif//dx, and making the change of variable f) = -2vln0(x, t), transforms the 
nonlinear Burgers equation into the linear heat equation. This provides the explicit solution of Eq.([T]i for any initial 
condition, and in the limit e — > 0^ a saddle-point method gives 



i/'(x, f) - min 
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'Paiq) 



It 



x-q(x,t) 

and v(x,t) , (6) 



where we introduced the Lagrangian coordinate q(x, t) defined by 

{x - qf 

i^o(q) -i — — is minimum at the point q-q(x,t). (7) 

The Eulerian locations x where there are two solutions q^ < q+ to the minimization problem dTji correspond to shocks 
(and all the matter initially between q^ and q+ is gathered at x), that is to clumps of particles of finite mass. The 
application q x{q, t) is usually called the Lagrangian map, and x q{x, t) the inverse Lagrangian map (which is 
discontinuous at shock locations). For the case of Brownian initial velocity that we consider in this paper, it is known 
that the set of regular Lagrangian points has a HausdorfF dimension of 1 /2 |25], whereas shock locations are dense in 
Eulerian space 1251 1241]. in the Brownian region. 

As is well known |6|], the minimization problem (|7]i has a nice geometrical solution. Indeed, let us consider the 
downwarcfl parabola fxA^) centered at x and of maximum c, i.e. of vertex (x, c), of equation 

(a - x)^ 

P.rAq)^-^-^+C- (8) 

Then, starting from below with a large negative value of c, such that the parabola is everywhere well below if/oiq) (this 

is possible thanks to the scaling i//Q(Aq) A^^^i//Q(q) for the integral i/tq of the Brownian motion, which shows that 
\i//o(q)\ only grows as q^^^ at large q), we increase c until the two curves touch one another Then, the abscissa of the 
point of contact is the Lagrangian coordinate q{x, t) and the potential is given by \l/{x, t) - c. 

This geometrical construction clearly shows that a key quantity is the conditional probability density, 
Kx,c{li' lAi ' ; 12, 4'!, ^'2)^ for the Markov process q i-» {iAo(^), VQ{q)], starting from {i/ri , vi ) at ^1 > 0, to end at V2) 



at q2 > q\ > 0, while staying above the parabolic barrier, if/Q{q) > PxAl), for q\ ^ q ^ qi- Following 111 111 (who 
studied the case of two-sided white-noise initial velocity) and [27] (who studied the case of two-sided Brownian initial 
velocity), we shall obtain the properties of the system from this propagator It was derived in [27J who obtained 

KxAqu^x,vuq2, ^2, V2) d(A2dv2 = G(t; ri, m,; r2, mz) drjdMz, (9) 

with 



T = r (62-61), r, =2y3 



Ui^2y(Vi + Qi-X). (10) 



Here we introduced the dimensionless coordinates (which we shall note by capital letters in this article) 

q X tv til/ tc I — 

2=-ir, X=— , y = -7, *I'=-7, C = ^, with y= V2Df, whence X^Q + V for regular points. (11) 

yl yl yl y4 y4 

For completeness, we give in AppendixlAlthe expression of the reduced propagator G and of two associated kernels. 

We can note that, thanks to the scale invariance of the Brownian motion, the scaled initial potential if^oi^q) has the 
same probability distribution as A^^^i//o(q), for any A > 0. Then, using the explicit solution Q we obtain the scaling 
laws 

i//(x, f) '= t'lj/ixlf-, 1), v{x, t) fv(x/f^ 1), q(x, f) '= t^qixlf-, 1), (12) 



'in the literature one often defines the velocity potential as v = -dyif/, which leads to upward parabolas. Here we prefer to define v = (9,1^ to 
follow 
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Figure 1: Geometrical interpretation of the initial conditions ipoiq") associated with the probability p^(0 < q' < 9, c)dc. The Brownian curve 
i//o(q") is everywhere above the parabola P^j. and goes below T^v.c+dc somewhere in the range < q" < q. On the left side, i//Q{q") = and it is 
also above the parabola P^ c- To obtain the cumulative probability, p^iQ < q' < q), we must then integrate over the height c of the parabola. 



where '^^ means that both sides have the same probability distribution. Indeed, the initial conditions on the left side, 
4io{q) - or ij/aiq) - +°°^ do not break these scalings. This implies that all our results can be written in terms of 
the dimensionless variables (fTTl i. as we shall check below. This would no longer hold if the Brownian domain were 
restricted to a finite interval [-L, L], since the size L would add a new length scale into the problem which can give 
rise to other scalings. In particular, at late times one would find a simple ballistic propagation |x| oc f to the left or right 
side in the empty case, when all high-velocity particles have already escaped from the Brownian domain. 

Finally, we may note that we defined the continuum limit as m — > at fixed density po. The same limit also 
describes the large time behavior of the system (i.e. f — > 00) at fixed reduced lengths and masses X - xjy^ and 
M - m/ipoj^), as in (fTTT i. In the first point of view, we consider the properties of the system at any finite time, arising 
from a distribution of infinitesimal particles, whereas in the second point of view we keep the discrete nature of the 
initial system but we look at the asymptotic late-time distribution over large scales and masses that grow as (so that 
discrete efi'ects become subdominant); see 11211 for more detailed analysis in the case of white-noise initial velocity. 



3. Case "F": expansion into an uniform medium at rest 

We first investigate the case "F", defined in (|4]l, where the left side, ^ < 0, is initially filled with particles at the 
same uniform density po with zero velocity. 

3.1. Eulerian velocity distribution Px{v) and probability of having been shocked ^^^ocked 

We first consider the one-point velocity distribution, Pxiv), at the Eulerian location x, as well as the distribution, 
Px(q), of the Lagrangian coordinate q(x, t) associated with the particle that is located at position x at time t. From 
Eq.® they are related by 

Pxiv) - t Pxiq) and q-x-vt, (13) 



since q(x, t) is well defined for any x except over a set of zero measure in Eulerian space associated with shocks ||24|] . 



As in 112711 . we first consider the cumulative probability, pxiO < q' < q), that the Lagrangian coordinate q'(x,t) is 
within the range Q < q' < q. From the geometrical construction (O, this is the integral over the parabola height c 
of the bivariate probability distribution, px(Q < q' < q, c)dc, that the first-contact point of the potential i/^o with the 
family of downward parabolas "Px.c, with c increasing from -00, occurs at an abscissa q' in the range Q < q' < q, with 
a parabola of height between c and c -i- dc. In terms of the propagator Kx^c introduced in Eq. (|9]l this probability density 
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reads as 



Px(0<q' <q,c)dc- lim I di//dvdil/-^dv+[K.c(0,0,0;q,ip,v) - Kxc+dc{0,0,0;q,il/,v)] Kxc{q,il^,v:,q+,il/+,v+). 

(14) 

Here we used the Markovian character of the process q i-> {iff, v), which allows us to factorize the probability p^(0 < 
q' < q, c)dc into two terms, which correspond to the probabilities that i) i//o stays above 'Px.c, but does not everywhere 
remain above 'Px.c+dc, over the range < q' < q, while reaching an arbitrary value {iff, v) at q, over which we will 
integrate, and ii) iffo stays above !Pv.e for q' > q. We show in Fig. [T]the geometrical interpretation of Eq. (fT4l i for a 
case with x > (we did not try to draw on the right side an actual Brownian curve if/oiq) which has no finite second 
derivative). 

The constraint associated with the left part of the potential i^fo at q < merely translates into an upper bound for 
the parabola height c. Thus, if x > 0, we must integrate over c up to the value c* - x^f(2t) where the parabola 'Px,c, 
runs through the origin {0, i/'o(0) = 0). Indeed, it is clear that all points on the negative axis {q < 0, iff^iq) = 0) are still 
located above 'Px.c. hence they cannot be the minimum associated with (|6]l. In fact, we can note that for any x > 
first contact always occurs before reaching c» because the initial potential has a zero derivative at q = (vo(0) - 0). 
This also means that no rarefaction interval opens at x = 0^ (nor at any other location, see On the other hand, 
for X < we must clearly integrate up to c = 0. Moreover, if the first contact is only reached at q - x for c = 0, it 
means that no particles from the right part have reached the position x yet. (For the case of two-sided Brownian initial 
velocity one would need to add a third factor of the form K^JQ, 0, 0; q-, iff-, v_) in Eq. (fT4l i to take into account the left 



part of if/o, see H27I1 .) 

Using the expressions given in Appendix A, as well as the results of appendices A and B of i27ll . we obtain from 
Eq.(fT4li for x > 0, after integration over c, the probability density 



r 

U — icx 



' zm 



which we expressed in terms of the dimensionless variables (fTTT i. Of course, the distribution vanishes over Q < 0, 
since particles from the left side cannot travel to the right side X >Q. 
In particular, at the origin X - this yields 

G > : Po(Q) = Q-''^ e-Q, whence P^m = {-Vy'l' with V < 0, (16) 

where we used V - -Q for X - Q. The probability vanishes for 2 < and V > 0, as particles cannot come from 
the left side. Thus we recover the results of Bertoin |3], who obtained Eqs.(fT5]l-(fT6]l from probabilistic tools. The 
distribution of the time increments of ^'(0, f), i e- of ^(0, 12) - q{0, ti), was obtained in fl. We can note that Eq.(fT6]l 
is identical to the large-velocity tail of the distribution obtained at x = for the case of two-sided Brownian initial 



conditions 112711 . Hence, for rare events the tail of the distribution does not strongly depend on the initial conditions on 
the opposite side of the origin. 

On the other hand, at large X we obtain for fixed velocity, V = X - Q, 

(,-v-ix 

X^+oo: PJV) ~ (17) 

Here we used the relationship Px(V) = Px(Q = X - V) between the probability distributions of the velocity V and 
of the Lagrangian coordinate Q, and the explicit expression ( fTsT l. Thus, as for the two-sided case, we recover at 
leading order the initial Gaussian distribution on large scales, here at x — » -i-oo. This is related to the "principle of 
permanence of large eddies" encountered in the hydrodynamical context [16], that holds for more general energy 
spectra, Eoik) cc k'\ with n <\. This states that regions of size x » L{t), where Lit) is the integral scale of turbulence 
(here L{t) ocy^ - 2Dfi, see the scalings (fTTTl). have not been strongly distorted by smaller scale motions yet (since the 
relative distance between particles has changed by an amount of order L(t) at time t). Thus, as checked in numerical 
simulations |1, 17], the stability of large-scale structures is not only a statistical property but actually holds on an 
individual basis, that is for each random realization of the velocity field. The properties of the velocity and Lagrangian 

6 



-0.5 -0.4 -0.3 -0.2 -0.1 O -2.5 -2 -1.5 -1 -0.5 

X = x/(2Dt2) log[-X] 

Figure 2: Left panel: The probability pj"*'='i that particles from the right side have already reached the position X on the left side. We show 
our results for the filled case "F" (solid line, Eq. j20) ) and the empty case "£"' (dashed line, Eq. )43t ). Right panel: Same as left panel but on a 
logarithmic scale. 

increments on the right part of the system, x > 0, were ah^eady obtained in and are also identical to those obtained 
far from the origin in [27] for two-sided initial conditions. In particular, it can be seen that the n-point distributions 
Px,,..,x,M\^ factorize as /?.vi(?i)p.v,,(?2,i)-Ai„,„-i(?«.«-i)' where we note Xij^i = Xi - Xj^i and = qi - the 
relative distances, for x\ < .. < x„ and qi < .. < q„. Thus, the increments of the inverse Lagrangian map, x q, are 
independent and have a simple distribution, which is given by the expression ( fTsT i without the factor i"'^"*. Then, over 
X > the properties of the density field p{x, f) are identical to those described in detail in [27] far from the origin, for 
the case of two-sided Brownian initial conditions. 

On the left side, x < 0, we must integrate Eq.(fT4li over c up to c = 0, as explained above. This yields 



J-/00 27n Jo -y/TT 

We can note that from Eqs.JTSIl. ( fTSl ), the tails at large Q and V read as 



-v2X + — 



over 2>0. (18) 



IX X 

Q^+oc: Px(Q)~— Q-^''^e-^, and, y ^ -oo ; PxiV) {-Vy^^'^e^, (19) 

^ r[i/4] ^ ' ' r[i/4] ^ ' ^ ^ 

which hold for both X > and X < Q. Thus, at any finite X the tail of the velocity distribution simply follows the 
exponential decay obtained at X = in ( fT6t . multiplied by a prefactor e^. 

The distribution ( fTSl l corresponds to realizations where some particles coming from the right side have already 
passed by position x. Thus, integrating Eq.lfTSb over Q >Q gives the probability, pj^^ed^ (.jj^j jjj^tjgj- coming from the 
right side has akeady passed by the position X <0. This yields 



Jo 



-v2X + \ 



(20) 



This could also be obtained directly by computing as in (fT4l l the probability that the curve i^o goes below the parabola 
!P.i_o at some point ^ > 0. Then, to the contribution (fTsT l we must add the contribution pno'-shocked^^^g _ 
pnot-,shocked _ ^ _ pshocked^ j-j^^j. corresponds to realizations where no particles from the right side have already passed 
by position x, so that the medium has remained at rest at x until time t. Finally, Eq.(l20li gives the asymptotic behaviors 



94/30 1/3 n4/3 \ c 

^ r[l/3] ^ 16 V -27tX 
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We can check that p^hocked - i at X = 0, while it shows an exponential tail at large negative X. Since all quantities can 
be expressed in terms of the scaled variables ( fTTT i. as we noticed from (T% and can be checked in the results above, 
this exponential decay can be obtained from simple scaling arguments. Thus, for particle ^ > to reach the Eulerian 
position X < Q at time t, we can expect its initial velocity to be of order vq ~ (x - q)/t. Since the initial velocity is 
Gaussian, of variance given by Eq.©, this corresponds to a probability of order g^'o/^^o?) - g-( t-9)-/(20f-i?) ^his is 
maximum at <7 = -x, which gives a weight _ g4x jjgjjcg recover an exponential over X. Of course, such 

a reasoning cannot give the numerical factor within the exponential. We show in Fig. |2]our results for the probability 
prfiocked^ given by Eq.(l20b. for the present case "F" where the left side is initially filled by particles at rest (solid line) 
(we also display for comparison the result associated with the empty case "£"', that we shall obtain below in Eq.(l43]l 
(dashed line)). We clearly see the fast decline with larger \X\ obtained in Eq. (l2T]) . 

3.2. Mean density profile and mean current 

From Eq.lfTsTi we also obtain the mean Lagrangian coordinate {Q{X)), and velocity {V(X)), at Eulerian location 
X > on the right side, as 

X>Q: <e(X)) = X+^, <y(X)) = -i, whence {q{x))^x+^-, {v{x))^~. (22) 

As expected, since particles gradually leak into the left side the mean velocity is negative, and particles that occupy 
the Eulerian position x come from increasingly far regions on the right as time increases. We also obtain the mass, 
m{< x), of excited particles (i.e. with initial Brownian velocity and which were initially located at ^ > 0) that are 
located to the left of the Eulerian position x by noting that it is given by m(< x) - poqix, t) since particles do not cross 
each other This yields from Eq.(l22l) 

x>Q: {m(< x)) - po{x + Dt^ /2), whence {p{x)} - po and {j{x)) - -poDt, with p — -r- and 

ox at 

(23) 

where we introduced the density p(x) and the current jix) at position x of excited particles. Therefore, we obtain a 
uniform mean flow from the right, with the mean current {j{x)) - -poDt, while the mean density remains equal to 
Po- Note that (jix)) - 2{p{x)}{v(x)), which implies that the fluctuations of the density and velocity are correlated. In 
fact, the velocity field v{x, t) is associated with Eulerian regular points, since shocks have a zero measure, but the flow 
of matter is associated with shocks since all the mass of excited particles is contained within shocks [24, 25, 5, 27ll . 
Therefore, it is not surprising to find that {jix)} + {p(x)){v{x)}, since these quantities probe different aspects of the 
dynamics. 

Thus, even though particles keep escaping into the left part x < 0, particles coming from the right semi-infinite 
axis keep replenishing the system and manage to maintain a constant mean density po over x > 0, through the mean 
uniform current -poDt that grows linearly with time. The linear growth with time of the mean velocity and current 
is due to the fact that at later times particles coming from more distant regions have been able to reach the boundary 
X = 0. Again, this exponent can be obtained from simple scaling arguments. Thus, at time t we can expect to see at 

the boundary, x = 0, particles coming from a distance q with an initial velocity of order vq q/t. Since the initial 

velocity scales as q^^^, see (O, this gives q ~ t^ and vq f, whence v f, assuming that matter flows through the 

boundary in a well-ordered fashion, so that the velocity of these particles has not been significantly damped by nearer 
lower-velocity particles, as the latter have already escaped into the left side. 

On the left side, X < 0, we obtain from Eqs.(fT8b-(l20l) the mean Lagrangian coordinate (QiX)}, and velocity 
<y(X)), as 

X<0: {V(X)) = -e^^ T ^3v-"^e-'^'" \(-X + y-^) Ai (-v2X + v'^) - v'^ Ai ' (-v2X + y-^)] , (24) 

Jo W V / V /J 

and {Q{X)) = X - {V{X)}. This gives for the mean velocity the asymptotic behaviors 



1 ?4/3t1/3 o [Zt^ 

X^O -. {V(X)} ~ -- + ^^(-X)"\ X^-oo: (y(X))~-- J— e^^/^. (25) 
^ ^ 4 r[l/3] ' ^ ^ 4 V ;r ^ ' 



-0.8 



-0.6 -0.4 -o.a 

X=x/(2Dt.2) 



-1.5 - 1 -0.5 

log[-X] 



Figure 3: Left panel: The mean velocity profile -2(v)/(Dt) on tlie left side, .v < (on the right side the mean velocity is constant and equal to 
-Dt/2). We show our results for the filled case "F" (solid line, Ea. )24t ) where the mean velocity is always meaningful. The dotted lines are the 
asymptotic behaviors 125 1 . Right panel: Same as left panel but on a logarithmic scale. 



We show our results for the mean velocity profile in Fig. [3] We can check that K^W)I shows a monotonic decrease 
for larger \X\ on the left side. As expected, we recover the same asymptotic exponential decay as for the probability 
^shocked obtained in Eq.dZB- 

The average (l24l i takes into account both the contributions (fTST l and f>™t-shocked^^g _ However, since only the 
part ( fTSl l contributes to the mass of excited particles located to the left of X, we no longer have (M(< X)) - {Q{X)), 
as was the case for X >Q 'm (|23] |. where we introduced the dimensionless mass M defined by 



M 



POT 



dim) d{M) 
2, whence <p) - — — - po and <;> 



dx 



dX 



d{m) 

'~dr 



.<M) 



(26) 



In the last two relations in ( l26b we used the property that {M{< X)) only depends on the reduced variable X. Integrating 
over the contribution (fTSl l gives 



X < : (M(< X)) = e 



2X 



r — 

Jo Vtt 



[v-^ Ai {-v2X + v"2) - v-^ Ai ' {-v2X + v"^)] 



(27) 



as well as 

jc < : <p(x)) = Po 
and 



„2X 



^3v-^''/2g-iv-'-v'x2 1(4^-3 _ f2X + - (4v"2 - 2vX) Ai ' i-v2X + v'^)] . 



X < : {j(x)) = -4pQDt e 



2X 



r dv 

— 3v 
Jo \n 



-3/2, 



[(v-3 _ 4v-3x + 6X^) Ai (-v2X + v~^) 



-(v-2 _ 4v-2x + 2vX^) Ai ' (-v2X + V"')] . 



(28) 



(29) 



We show in Figs. |4] |5] our results for the mean density and current. Then, we can check that the mean density and 
current are continuous at the boundary x -Q, and we obtain the asymptotic behaviors 



1 2^/332/3 

X^Q- : <M«X))~ -+X-^J^(-X)5/^ <p(x))~po 



95/30-1/3 



{j{x)} ~ -poDt 



714/30-1/3 



(30) 



-0.6 -0.4 -0.2 

X=x/(2Dt.2) 
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Figure 4: Left panel: The mean density profile {p)/po on tlie left side, x < (on the right side the mean density is constant and equal to po). We 
show our results for the filled case "f " (solid hne, Eq.(28j) and the empty case "£"' (dashed line, Eq. )46t ). Right panel: Same as left panel but on 
a logarithmic scale. 




Figure 5: Left panel: The mean cuirent profile -{]) lipoDt) on the left side, .v < (on the right side the mean current is constant and equal to 
-poDt). We show our results for the filled case "F" (solid line, Eq. )29t ) and the empty case (dashed line. Eg. 1471 ). Right panel: Same as left 
panel but on a logarithmic scale. 
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From Eq.(l30b. we can note that the amphtude of the mean density and current first increases to the left of the 
boundary x - Q (contrary to the mean velocity which showed a monotonic decrease into the left side). This is due 
to the dragging effect of the matter that was initially at rest on the left side, which slows down the leftmost cluster 
associated with particles coming from the right side. Then, this deterministic friction leads to a transitory growth for 
the mean density and current to the left of the boundary x - 0. Figures |4] and |5] show that although this feature can be 
clearly seen for the mean density (although it remains modest, of order 5%), it is almost invisible for the current, in 
agreement with the higher power (-X)^^^ instead of (-X)^^^ in Eq. dSOl) . which leads to a suppression by a factor \X\ 
at small X. As we shall check in section l4!2l below, this feature disappears when the left side is initially empty. The 
same behavior was obtained in (]3\ for the case of white-noise initial velocity. At large distance from the boundary 
we obtain the exponential decays 

(31) 

Again, the characteristic length scale is the reduced variable X of (fTTT i. and the exponential decay is the same as the 
one obtained in Eq.dllll for the shock probability pj°=ked 

3.3. Lagrangian displacement field 

We now consider the dynamics from a Lagrangian point of view. Thus, labelling the particles by their initial 
position q at time f = 0, we follow their trajectory x(q, t). Since particles do not cross each other, the probability, 
Pq{x' > x), for the particle q to be located to the right of the Eulerian position x at time f, is equal to the probability, 
Pxiq' ^ q), for the Eulerian location x to be "occupied" by particles that were initially to the left of particle q. In terms 
of dimensionless variables, this gives for right-side particles, q>Q, 

2m 5-1 

where the integration contour runs to the right of the pole, i.e. ^R(i) > 1, and 

e > : P,iX' > X) . 1 - r ^'^r i;3v-3/^.-|.^''^--.,.^. 

J^iac 2m 1 - s Jo 

where the integration contour obeys < %{s) < 1 . Note that for X < we must take into account both the probability, 
1 _p^°';ked^ jjj^j- particles from the right side have reached X yet, and the probability, Px{0 < Q' < Q), that particles 
Q' from the right side, with < Q' < Q, have already passed by point X. Then, probability densities are obtained 
from Eqs.(l32b-(l33]) by differentiating with respect to X. 



-v2X + — 



over X<Q, (33) 



3.4. Leftmost cluster 

We can identify the Eulerian position, ji:(0, t), of the particle that was initially located at the origin (q = 0), as the 
position of the leftmost cluster (or "leader") formed by excited particles that have escaped into the left side. In the 
present case, where these particles spread into a medium of uniform density that was initially at rest, this cluster also 
contains a mass pol-^l of particles that were located in the interval [x, 0]. It acts as a snow-plough while the conditional 
properties of the system to its right are no longer sensitive to the initial conditions on the left side. From Eq.(l32b we 
can see that Po(X) vanishes for X >0, which means that at any time f > the particle ^ = has almost surely already 
passed to the left side and the leftmost cluster has formed with a finite mass, in agreement with the results of H. On 
the other hand, for X <Q the probability density PidX) of the leftmost cluster position reads as 



X<0: Pie (X) = Apshocked ^ _d_ ^2X T iL3^-3/2^-|v--v'X=^j 

' dX ^ dX Jo 



-v2X + 4r 



(34) 



Here we used the obvious property PidX' > X) - pnot-^hocked _ j _ pshocked^ which states that the leftmost cluster 
is located to the right of point X if, and only if, no particles from the right side have reached this point yet. We can 
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Figure 6: Left panel: The probability distribution P\j^_(X) of the position X of the leftmost cluster. We show our results for the filled case "F" (solid 
line, Eq. )34t ) and the empty case (dashed line, Eg. 1501 ). Right panel: Same as left panel but on a logarithmic scale. 



check that the result ( l34b is identical to the one that would be obtained from Eq.(l33Tl. Using Eq. dSTT i we obtain the 
asymptotic behaviors 



PiAX) 



24/3 3-2/3 



-i-xy 



-2/3 



MX/9 



(35) 



r[l/3] " -....V". 32/3 V _2;rX 

We show in Fig.|6]our result (|34] | for P\.c.(X). From Eq.(l34b. we obtain after an integration by parts the mean position 
of the leftmost cluster as 



{xiAt)} = <X,.e.)2Df2 with <X,,.) 



r" ,^ dv 

dXe^^ — 3v 
J-00 Jo w 



-3/2, 



"Ai 



-v2X + 



1 



^ -0.06. (36) 



Thus, the distance |jci.c.(f)l from the origin scales with time as f. The leftmost cluster has a super-balUstic motion 
because it is constantly overtaken by higher- velocity particles coming farther away from the right side which increase 
its momentum. Since the growth of the mean {xidt)} is set by the scaling variables (fTTl i it can be understood from 

simple arguments, as for the mean current obtained in ( l23T l. There, we have seen that the scaling v'(0, t) 1 could be 

explained by the time needed for particles at distance q ~ p- to reach the origin. Then, if the position of the leftmost 

cluster is set by the latest particles that escaped into the left side we expect xi.c. ~ vt which agrees with ( |36] |. 

Next, to derive the distribution, p\.c.{m), of the excited mass of this leftmost cluster, we first consider the bivariate 
distribution, pic.(x, < q' < q)Ax, that this aggregate is located at a position in the range [x, x + dx], with a right 
Lagrangian coordinate q' that is smaller than q. This corresponds to a mass m' - p()q' that is smaller than m - poq. In 
a fashion similar to ( fT4] i. we can write this quantity in terms of the Brownian propagators K^ c as 



pic.(x, < q' < q)dx ■ 



Hm I diAdvdiA+dv+ [^^1,0(0, 0, 0; q, 41, v) - K^+i^fli^, 0, 0; q, iff, v)] K,fi(q, if/, v; q+, if/+, v+). 

(37) 

This equation states that the initial velocity potential if/o obeys the following two constraints, i) it stays above !P , but 
goes below P_t+dx,Q, over the range < q' < q, and ii) it remains above !Pj,o for q' > q. Using the expressions given 
in Appendix A, as well as the results of appendices A and B of [27.1 . we obtain from Eq.(l37]i. in terms of the reduced 
variables X and Q, 

PiAX,0<Q'<Q) = 2e2x I 3^- 

J_,„ 2m s - 1 Jo 



[^/^-v^x)Ai 





- vAi' 






-v2X + 4 


-v2X + 4 


} 











with %(s) > 1. (38) 
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Figure 7: Left panel: The probability distribution Pi.c.(M) of the mass M of the leftmost cluster We show our results for the filled case "F" (solid 
line, Eq. )39t ) and the empty case (dashed line, Eg. 1551 ). Right panel: Same as left panel but on a logarithmic scale. 



We can check that the Umit Q — > +oo, which is governed by the pole at i = 1 in Eq.dSSll. yields back the probability 
density Pi.c.(^) given in Eq.(l34ll. Integrating over X, using an integration by parts and differentiating with respect to 
Q gives the mass distribution 



r[l/4] X,.^ Ini ' Jo Vi 



-v2X + 4 



(39) 



where the reduced mass is given by M = m/ipoy ) = Q, as in ( |26] |. This yields the asymptotic behaviors 



M ^ 0+ : Pi.c.(M) ~ f^J^ M-^'*, M ^ +00 : PiAM) ~ aM'^'^e'^ (40) 

with a = [2i°/^3-i^''/^2fi(f,5;f;|) + 223/33-17/3 2Fi(|,f;f;i)]/V^ ^ 0-511- We show in Fig.Elourresuh (O 
for the probability distribution of the mass of the leftmost cluster. We can note that the high-mass tail has the same 
forrn, M"3/2e-'*^, as the mass function of shocks on the right side x > (but with a slightly smaller normalization), 
see 0311271]. From Eq.(l39]l. integrating over M and s, we recognize the integral in Eq.(l36ll and the mean excited mass 
of the leftmost cluster reads as 

<mi.e.(f)> = <Mi,.)po2Z)f2, with <Mu.) = ^ + <X,,.) - 0.19. (41) 

Thus, the mass of the leftmost cluster grows with time as f in the mean, as it keeps being overtaken by new particles 
that arise from more distant regions {q ~ fi) on the right side. We can note that the particles that were initially at rest 
on the left side, and slow down its propagation, do not change the scaling laws associated with this cluster, as shown 

by Eq.(l36ll. Indeed, since xic f, the leftmost cluster has captured a mass rrf'^°^^'"^ ~ of still particles that were 

initially at rest on the left side. On the other hand, it has just been overtaken by a mass m ~ P' of excited particles 

with q ~ f- and v 1. Therefore, since the mass nf'^°'^^^^ scales as m, the total momentum of the leftmost cluster 

still scales as the momentum of the excited particles. As a consequence, we shall obtain in section l431 below the same 
scalings for the case "£"' where the Brownian particles expand into empty space. 

4. Case "£■": expansion into an empty medium 

We now consider the case "£"', defined in (|5]l, where the left side, <7 < 0, is initially empty. We also compare our 
results with the previous case "F". 
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4.1. Eulerian velocity distribution and shock probability 

From the fact that the initial velocity potential is infinite on the left side, i/'o(^) = +°° for ^ < 0, we can see that 
the computation (fT4l l-(fT5]l. performed in section lTTI for the case of a filled left side at rest, remains valid in the present 
case. Moreover, contrary to the previous case, this same computation also holds on the left side, x < 0, since the 
highest height c, of the parabolas !Pv,c is still set by the origin [q = Q,if/o = 0). Therefore, for all X we obtain the 
probability distribution of the Lagrangian coordinate Q as: 



ioo Jo y7T 



Ai 



-v2X + — 



over g > 0. 



(42) 



Of course, the Lagrangian coordinate Q is always positive, since the left side is initially empty. Note that the only 
difference with Eq.lfTS]) is the absence of the factor e"'' that arose from the additional constraint PxAo) - o^^r 
< in the filled left-side case "F" for x < 0. For X > 0, we can integrate over v, using the results of Appendix B 
in ll27[l . which gives back Eq.lfTsll. Thus, all properties of the system on the right side, X > 0, are identical to those 
obtained in section|3]for case "F". Moreover, at large Q and V, we recover the same asymptotic behaviors ( fT9] ). 
The shock probability, p^hocked^ jj^^j. g^cited particles have already reached the position X < on the left side, now 



reads as 



X < : P 



ishocked 



= e 



2X 



dv 



3v 



-3/2, 



Ai 



-v2X + — 



(43) 



Again, it only diff'ers from Eq.(l20ll by the absence of the factor e^"^^^ . Since the Airy function Ai(3;) is positive over 
y > this shows that p^°'=ked^ 2& given by Eq.(l43]l. is always larger than the value (l20l i obtained in section [TTI for the 
case "F". This is due to the fact that in that previous case "F", particles that escape from the right side are slowed 
down as they travel to the left by the matter that was initially at rest on the left side. (More precisely, this matter slows 
down the leftmost cluster, which also slows down particles that overtake it and aggregate to it.) This clearly implies 
that, at any location X < 0, is smaller for the filled left-side case "F", as can be checked in Fig. |2] where we 

compare the result (l43T l (dashed line) with the previous result (|20] | (solid fine). From Eq.(l43]l we obtain the asymptotic 
behaviors 

-1A-X 1 

1 - ^/ , X ^ -00 : p^^hocked ^ __ g6X_ (44) 



x^Q- 



[jshocked 



We can check that they are larger than the results (ISTT i obtained for case "F", in agreement with the discussion above. 
We can note that the exponent 1/3 has been changed to 1/2 for the limit X — > , whereas the large-X decay still 
has the form of an exponential multiplied by an inverse square-root, but with a different numerical factor in the 
exponential. 

4.2. Mean density profile and mean current 

As discussed in section HTTl the properties of the system, whence the mean density and current, are identical over 
x > to those obtained for the previous filled left-side case in section|3] In particular, we recover the results (l22l i-(l23Tl. 
Therefore, we now focus on the left side x < 0. For the present case, where there are no particles at x < until some 
particles from the right side have travelled down to position x, the mean Lagrangian coordinate {q{x)) and the mean 
velocity (v(x)) are not meaningful. However, we can still define the mean mass (;«(< x)) of excited particles located 
to the left of the Eulerian position x, as well as the mean density and mean current (all these quantities being equal to 
zero until particles from the right semi-infinite axis have reached position x). Then, we obtain from Eq.(l42li 



X < : <M(< X)) = e 



Jo \n 



-3/2, 



[v^^ Ai (-v2X + v"2) - Ai ' (-v2X + v"^)] , 



X < : <p(x)) 
and again for x < 0, 



„2X 



r iL3v-3/2g-§v- [(4^-3 _ 4^) /_^2X + - 4y-2 Ai ' {-v2X + y-^)] . 
Jo w 



(45) 



(46) 



{j(x)) = -4poDte- 



2X 



Jo \7T 



[(v-3 _ 4y-3x + AX^) Ai {-y2X + y-^) - (v^^ 



■Ay-^X)A\' {-v2X + v-^)\. 

(47) 
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This yields close to the boundary x - the asymptotic behaviors 



1 2''/23'/2 
(M(< X))~-+X+ ——(-Xf\ 
4 5 

2ll/2o3/2 

i-^—^(-xf' 

5 yn 



{p(x)} ~ po 



27/23I/2 



-(-X) 



3/2 



{j(x)} ~ -poDt 



(48) 



Thus, we can see that the transitory increase obtained for the amplitude of the mean density and current in the case of 
the filled left side, see Eq.dSOll. that was associated with the dragging effect due to matter that was initially at rest on 
the left side, is no longer present when the Brownian particles expand into an empty medium. At large distance we 
obtain the same exponential decay as for the shock probability pj°<:ked^ 



-co 



<M(< X)) 



-2X 



fix 



{p(x)) ~ pq2 



-6X 



fix 



{j(x))--pQDtS^I-{-Xr"e' 



3/2 ^6X 



(49) 



Again, we can check that, far from x = 0, the amplitude of the mean mass, density and current are larger than for the 
filled left-side case given in dSTT i. This agrees with the larger value of p^°<:ked discussed above. We can check these 
properties in Figs.|4]|5] where we compare the mean density and current obtained for both cases. 

We can note that the propagation into empty medium for the Brownian case studied here shows significantly 
different properties from the case of white-noise initial velocity studied in 1 13]. Indeed, in the latter case, the typical 
distances and masses only scale as fi^^, in a fashion similar to the scaling associated with the scaling laws (fT2l l 
that express the scale invariance of the initial velocity field. However, the mean mass, (M(< X)}, located to the left 
of X, only decays as the inverse power-law 1/X^ (with now X ~ xjf-l^ and M ~ m/t^^^), instead of the exponential 
falloff' (|49T l. Then, one still obtains a finite mass distribution in the ballistic limit f — » 00 for p^fyn, t) at fixed ^ and m. 
This corresponds to free-moving forerunners that carry a finite mass in front of the typical profile that extends on the 
smaller scale ~ fi^^ . For the case when particles spread into a filled medium at rest, this distribution vanishes [Tsl]. 
In the present Brownian case, the appearance of such a new nontrivial scaling, specific to the expansion into empty 
space, is no longer possible, as shown by the exponential decay (|49] | and the results derived above. Indeed, the front 
now shows a faster than ballistic propagation in both cases "F" and "£"' - as as given by the scalings ( fTSl l - which 
is governed by the latest high-velocity particles coming from increasingly far regions on the right (see the discussion 
below Eq.dTTTl). Then, this super-ballistic propagation leaves no room for new scalings and the dynamics remains very 
similar for both cases "F" and "£"' (see also 1119(1 for a comparison of white-noise and Brownian initial conditions). 



4.3. Leftmost cluster 

Within a Lagrangian point of view, the probability distributions of the displacement field are again given by Eq.(l32li 
overX > 0, and by Eq.(l33]l over X < but without the factor e in agreement with the previous discussion. Then, 
we focus here on the leftmost cluster, which is also associated with the particle ^ = 0. From Eq.(l43Tl we now obtain 



X<0: Plc(X)= pshocked 

dX 



which leads to the asymptotic behaviors 
x^O--. PiAX) 

while the mean position is 



dX Jo 



-nX 



-"^e-i^-'Ai 



-v2X + — 



X 



Plc.iX) 



-nX 



„6X 



(50) 



(51) 



X r°° dv 
dXe^^ —3v-^'^e- 
00 Jo 



"Ai 



-v2X + 4 



^ -0.08. 



(52) 



As discussed in section [341 we recover the same scaling law as for the case where Brownian particles expand into a 
filled medium at rest, see Eq. (l36t . However, the mean of the reduced variable Xi^. has a slightly larger absolute value, 
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since it is easier for the leftmost cluster to travel to the far left as it is no longer slowed down by the particles that were 
initially at rest in the filled case "F". This also leads to the smaller low-X tail and to the larger high-X tail that can be 
seen in Fig. |6] and by the comparison of Eq.dSTTl with Eq.(l35Tl. Note however that for the most part both distributions 
P\.c.(X) are very close to each other 

The bivariate distribution p\.c.{x, < q' < q)dx is now given by 

PiAx, 0<q' < q)dx = ^ lim^ J" d(Advd(A+dv+ [K^^^2/(2r)(0, 0, 0; q, v) - /r.v+dx,x2/(20+(x/0dx(0, 0, 0; q, i//, v)] 

X ^x,xV(2/)(?, lA, v, q+, v+). (53) 

Indeed, since tfroiq) = +oo over < 0, we no longer have to consider parabolas V^.c with c - 0, which are tangent to 
the horizontal axis i^o = 0, as we did in Eq.dSTli for the case "F". In the present case "£"', we must consider parabolas 
which run through the origin, !Px,c(0) = 0, whence c = x^/(2t). Using the expressions given in Appendix A, as well as 
the results of appendices A and B of we obtain from Eq.( 



27n s - 1 Jo I 



-v2X + — 



- vAi' 



-v2X + — 



(54) 

where the integration contour obeys %(s) > 1. Again, taking the limit Q — > +oo we can check that we recover the 
distribution (l50l l. Then, integrating overX yields the mass distribution 



yV/-3/4 



r[l/4] 2m 
This gives the asymptotic behaviors 



J-,oo 27n J_oo Jo yj^ 



-3/2 -f.^ 



Ai 



-v2X + 



(55) 



M -» 0+ 



Pic.(M) 



V372 

r[i/4] 



and the mean mass 



1 



(m,,.(f)> = {Mi,,,)po2Dt\ with (Mi.e.) = - + (Xi,,,) - 0.17, 



(56) 



(57) 



where we recognize the integral ( |52] |. in a fashion similar to (1411 1. Thus, as for the position xidt), we recover the 
same scalings for the mass mi c.(f) as for the case "F" of the expansion into a filled medium at rest. The mean reduced 
mass (Ml c ) is slightly smaller than for the case "F". Indeed, since the leftmost cluster is no longer slowed down by 
particles that were initially at rest on the left side, it moves somewhat farther into the left side, as seen in Eq.(l52li. 
which implies that fewer particles from the right side have been able to overtake it. This leads to a smaller mass of 
excited particles that have been able to aggregate into this cluster. This now implies a larger low-M tail and a smaller 
high-M tail, as compared with the case "F". These properties can be checked in Fig. [T] and by the comparison of 
Eq.(l56t with Eq.(l40ll. However, as for the position Xi^., both distributions PidM) remain for the most part very close 
to each other 



5. Conclusion 

In this article we have studied the one-dimensional ballistic aggregation process, in the continuum limit and for the 
case where the initial velocity on the right semi-infinite axis is a Brownian motion. The left side is either at rest, with 
the same uniform initial density, or empty. Then, focussing on the out-of-equilibrium propagation of particles towards 
the left of the system, we noticed that in both cases the mean density remains constant on the right side whereas a 
mean current towards the left develops and grows linearly with time. Thus, particles coming from increasingly far 
regions on the right replenish the system, as seen at any finite distance x > to the right of the origin, and balance 
the mean loss of matter associated with particles that have escaped into the left semi-infinite axis. Moreover, the 
properties of velocity increments, of the density field and of shocks are the same in both cases and are also identical 
to those obtained asymptotically far from the origin for the case of two-sided Brownian initial conditions. 
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We find that on the far left, for both cases all quantities (e.g. mean density and current) show an exponential decay, 
with a slightly faster falloff for the case "F" where the particles expand into a filled medium at rest. Indeed, in this 
latter case, the initially still particles that were located on the left side slow down the propagation towards the left of 
the leftmost cluster, built by the boundary particle, ^ = 0, and by all particles ^ > that have already overtaken it, as 
they aggregate to it and decrease its total momentum. This dragging effect also leads to a transitory increase to the left 
of the boundary x = of the mean density and current, which is not present when particles expand into empty space. 
This also leads to a leftmost cluster which is statistically closer to the boundary x - Q and more massive (counting 
only the particles that came from the right part) as it is easier for particles from the right side to overtake it. 

An interesting feature of this one-dimensional process is that it provides a nontrivial inhomogeneous non-equilibrium 
system where many quantities can be obtained explicitly, as seen for instance in the calculations presented here. The 
same methods could also be applied to different time statistics, but we leave such studies for future works. Then, 
the simple system described here may also be used as a benchmark to test approximation schemes devised for more 
difficult cases where it is not possible to derive exact results. 

A. Brownian propagators 

The reduced propagator G introduced in Eq.® is most easily written in terms of its Laplace transform, G, defined 
by a 

/ 

Jo 



G(s;r\,ui;r2,U2) - I dre ^^G{T;r\,u\;r2,U2), and G-Gq-G\ with 
io 
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The term Go actually corresponds to unconstrained Brownian trajectories, hence it is also given by 10] 
r l-r- ^ ,, ■ ^ ,, \ — ^ „-4-(''2-n-"iT)- + 4-(''2-n-»iT)(ti2-tii)-jr(»2-Mi)- 

(jq(T, ri,U\,r2,U2) - ^^^2 ^ r- r- r _ 



[-0(-v)0(n - rz) + e(v)0(r2 - ri)] , (59) 
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From the propagator G it is convenient to derive the kernel H^o , associated with Brownian particles that remain forever 
above the parabola P^^, 



lim A\l/2Av2K^_c{qu^\,vx;q2,\p2,V2) - e "''^H^iruui). 



(62) 



We also consider the propagators, A and E, associated with Brownian particles that come within a small vertical 
distance 6c, or horizontal distance 6x, from the parabolic absorbing barrier: 

lim ^[/r.i-,c(?i, lAi. ^i; ?2, V2)-K_,,c+sciq\,i/^i,vi;q2, ^2, V2)] di/rzdvi = 2- e-^/r-+(«2-«i)/r mi; r2, U2) dr2du2, 

Sc^o oc y 

(63) 

and 



lim — {K^,c{q\, lAi ' ; ?2, 1A2, V2) - K^+6xA<i\ , lAi . ; ?2, v-a)] diA2dv2 = 

5a— >o ox 



-1 ^-T/r +(«2-«i)/r 



E(t; ri, uu r2, U2; qi, q2, x) dradMa. 



(64) 



Using Eqs.(l59]l-(l60]l, one obtains the expressions 
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